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OF  SEVERAL  EXPONENTIAL  LOCATION  PARAMETERS 

Hubert  0.  Chen  * 

Department  of  Statistics  and  Computer  Science 
University  of  Georgia 

ABSTRACT 

Let  there  be  k  independent  exponential  populations  with  unknown 
a;,d  possibly  unequal  location  parameters,  and  with  a  common  but 
uiiKnown  scale  parameter.  A  complete  random  sample  of  common  size 
n  is  drawn  from  each  population.  Based  on  these  k  samples,  this 
paper  is  to  propose  a  new  range  statistic  called  C-statistic,  for 
comparing  these  exponential  location  parameters  by  a  set  of 
simultaneous  confidence  intervals.  The  upper  percentage  points 
of  the  C-statistic  are  given,  and  the  power  function  of  the 
statistic  under  the  alternative  hypothesis  is  also  studied.  Other 
applications  in  multiple  range  test  and  in  ranking  and  selection 
problems  could  be  expected. 
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1.  INTRODUCTION 

In  life  testing,  there  are  many  occasions  when  the  exponential  distributions  are 
appropriate  statistical  models  for  the  random  variables  under  study.  The  problem  of 
testing  the  hypothesis  concerning  the  parameters  of  these  distributions  has  been  ex¬ 
tensively  studied  in  the  past  four  decades.  Paulson  (1941),  Epstein  and  Tsao  (1953), 
Kumar  and  Patel  (1971),  Weinman,  et  al.  (1973),  and  others  considered  the  inference 
Gri  two  exponential  location  parameters  based  on  likelihood  ratio  or  its  equivalence. 
For  testing  the  equality  of  several  exponential  parameters  Sukhatme  (1936)  considered 
likelihood  ratio  tests;  Hogg  and  Tanis  (1963)  proposed  an  iterated  testing  procedure 
by  which  earlier  conclusion  may  be  drawn.  These  procedures  are  long  well  known  and 
only  applicable  for  testing  a  hypothesis.  However,  the  multiple  comparisons  between 
these  parameters  of  interest  by  either  pairwise  differences  or  a  linear  contrast  were 
not  reported.  In  this  paper,  we  propose  a  new  range  statistic,  C,  which  is  similar 
to  the  Studentized  range  statistic,  for  comparing  several  exponential  location  para¬ 
meters  by  their  differences  or  a  linear  contrast. 

The  statistic  C  proposed  in  this  paper  has  not  been  seen  in  the  literature  nor 
in  the  books  by  Gumbel  (1958),  Sarhan  and  Greenberg  (1962)  and  David  (1981).  The 
distribution  of  C  and  its  moments  are  therefore  obtained  in  Section  2.  Applications 
of  the  C  statistic  to  hypothesis  testing  and  multiple  comparisons  concerning  these 
location  parameters  are  given  in  Section  3.  The  power  function  of  the  test  is  con¬ 
sidered  in  Section  4.  And  the  table  of  upper  percentage  points  is  also  provided. 


■sajaaasaaaSteaagiis 
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2.  DISTRIBUTION  OF  C 

Let  X^,  ....  be  independent  and  identically  distributed  random  variables 
(l.l.d.r.v.'s)  with  the  standard  exponential  probability  density  function  (p.d.f. }> 
f(x)  =  exp(-x),  x  >  0,  zero  elsewhere.  Let  R  =  max(X.. )  -  min(X^)  be  the  range  of 
the  X's.  Then,  it  is  easy  to  see  that  the  distribution  of  R  has  the  p.d.f.  of  the 
form 

g(r)  =  ( k~ 1 )exp(-r ) { 1  -  exp(-r)}k~2,  r  >  0,  (2.1) 

zero  elsewhere.  Let  V  be  a  chi-square  r.v.  with  2v  degrees  of  freedom  (d.f.),  in¬ 
dependent  of  the  X's.  Define  the  range  statistic  C  to  be 


C  V/2v  * 

Let  Y=R.  Then  the  joint  p.d.f.  of  C  and  Y,g(c,y),can  be  easily  found  as 

g(c.y)  -  (k-l)vvyVyU+v/i:)(l-e->')k-2/(r(v)cv*1). 

c  >  0,  y  >  0, 


(2.2) 

(2.3) 


zero  elsewhere.  The  marginal  p.d.f.  of  C  can  then  be  obtained  by  using  the  binomial 
expansion,  collecting  the  exponential  terms  and  integrating  out  y  in  g(c,y),  we  have 
the  p.d.f.  of  C,  h(c),  where 

k-2 

h(c)  =  (k-1)  E  ('l)i(k:2)[l+c(i+l)/v]"(v+1),  c  >  0,  (2.4) 

i=0  1 

zero  elsewhere. 

We  note  that,  if  k=2,  the  C  has  an  F  distribution  with  2  and  2v  d.f.  and  it  has 
an  L-shaped  density.  When  k>2,  the  distribution  of  C  has  an  unimodel  occurring  at 
c  =  cQ  >  0  and  the  shape  of  the  density  is  skewed  to  the  right. 
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The  upper  ath  percentage  point  c  a  c,a  is  obtained  according  to  the  equation 

K  |  v 

k-2 

C  h(c)dc  =  s  (-l)i(i[:}){l+c(i+l)/v}‘v  -  a.  (2.5) 

c  i=0  11 

Since  the  range  statistic  will  be  used  in  certain  applied  problems  associated 
with  samples  of  size  n,  we  let  v  =  k(n-l)  and  then  compute  the  table  of  the  percentage 
point  c  for  k  =  2(1)10(2)20,25,30,40,50,  n  =  2(1)10, 12, lb, 20, 30, 60, 120,»;  a  =  .10, .05, .01 
which  is  given  in  Table  1. 

The  mth  moment  of  C  about  origin  may  be  obtained  by  the  changes  of  variable 
y  -  l/(l+c(i+l)/v).  Thus,  we  have 

E(Cm)  =  r  cmh(c)dc 

k-2  .  .  r,  m+1  m  .  , .  . 

=  (k-1)  I  (-l)1(kT2)(^T)  /'  l  (f)(-l)Jyv*m“1+J4y 

i=0  1  11  0  j=0  J 

k-2  L  9  nr*- 1  m  . , . 

=  (k-1)  £  (  i  )(7TT)  £  (-1)1  ^(^/(v-m+j) 

i=0  1  11  j=0  J 

provided  v >  m. 

It  is  easy  to  show  that,  as  n  or  v  goes  to  infinity,  the  distribution  of  C  in 
(2.4)  converges  to  g(r).  Thus,  the  critical  value  of  c  can  be  obtained  by  the  equation 

c  -  -tna-a-v)1^-1*)  .  (2.6) 

For  large  n,  the  approximation  is  fairly  good  as  compared  with  the  tabled  values. 


-  - tggjj 
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3.  MULTIPLE  COMPARISONS 

Let  tt j ,  tt2,  ....  denote  k(>2)  populations  so  that  n  independent  observations 
Xil*  X-j2*  •••»  X^n  taken  from  population  are  exponentially  distributed  with  the 
probability  density  function  (p.d.f.) 

f(x;oti,e)  =  (l/0)exp{-(x-ai  )/0} ,  ot^  <  x  <  0  >  0,  (3.1) 

zero  elsewhere,  where  ot^  is  an  unknown  location  parameter  and  6  is  a  common  but  un¬ 
known  scale  parameter  for  i  »  1,  2 . k.  It  is  known  that  the  location  parameter 

is  also  referred  to  as  the  guaranteed  life  span  and  the  scale  parameter  the  standard 
deviation. 

Let  Y..  denote  the  first  order  statistic  of  the  sample  of  size  n  from  population 
n.j .  Let  Y^  <  ...  <  Y^  denote  the  ordered  values  of  the  k  first  order  statistics 
Yj,  ...,  Yk  and  let  the  well  known  minimum  variance  unbiased  estimator  of  e  be  given 
by 

k  n 

e  -  E  E  (x,,  -  Y,)/(k(n-l)). 
i=l  j=l  13  1 

Define  the  range  statistic  C  to  be 


c  =  n(Y[k]  "  Y[l]^®- 


(3.2) 


Under  the  null  hypothesis  HQ  :  oij  =  . . .  =  ak,  C  has  the  same  distribution  as  of  (2.4). 

A 

This  is  because  C  is  the  ratio  of  R  =  n(Y^  -  Y^)/etoV  =  2v0/0,  where  R  and  V  are 
independent. 

The  C  range  statistic  provides  a  quick  test  for  Ho;  one  rejects  HQ  at  a  level  of 
significance  if  the  computed  value  of  C  is  larger  than  c.a  where  P(C>c.a  |  H  )  -  a. 

K  t  V  K  y  V  O 
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Since  the  C  statistic  is  based  on  the  range,  it  will  not  be  as  powerful  as  the  F-test 
considered  by  Sukhatme  (1936)  when  a  test  of  hypothesis  is  concerned.  However,  if 
the  multiple  comparison  among  a's  is  interested,  the  Sukhatme's  test  is  not  applicable. 
Like  Tukey's  Studentized  range  statistic  for  comparing  normal  means,  the  major  role 
of  the  C  range  statistic  will  rest  on  its  extensive  use  in  pairwise  multiple  comparisons 
and  simultaneous  inferences  on  the  location  parameters. 

Similar  to  Tukey's  pairwise  multiple  comparison  for  normal  means,  we  can  construct 
a  set  of  simultaneous  confidence  intervals  for  all  of  the  differences  a.  -  a ■  by  the 

*  J 

following  probability  statements: 

P<HW  -  -  ck,k(n-l)  ’  n  *  “  1 . k-  ‘  1  ‘  “• 

This  is  because  the  inequalities 


(W  -  <°rV  ,  e 

8/n 

hold  true  for  all  i  f  j  if  and  only  if  the  inequality 


(3.3) 


§  max  {|(Yi-ai)  -  (Yj-Cj ) I > 

_ _ _ 

e/e 


<  c 


(3.4) 


holds  true.  The  numerator  in  the  left  member  of  the  inequality  (3.4)  is  the  range  of 
k  independent  standard  exponential  r.v. 's  and  the  denominator  is  the  chi-square  r.v. 
(with  2vd.f. )  divided  by  2v.  Their  ratio  defines  C  as  of  (2.2).  Therefore,  a  set  of 
(l-a)100%  simultaneous  confidence  intervals  for  the  difference  ai  -  is  given  by 


ai  -  aj  e  (Yi  -  V  1  ck?k(n-l)  ‘  5/n 


(3.5) 


for  all  i,j  =  1,  2,  ....  k,  i/j.  In  the  case  where  the  sample  sizes  are  not  all 
equal,  we  suggest  to  replace  1/n  in  (3.5)  by  (1/n.j  +  l/nj}/2  and  k(n-l)  by 
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k 

v  =  Z  (n.-l),  in  which  a  conservative  set  of  confidence  intervals  is  conjectured. 
1  1 


Furthermore,  a  generalization  of  the  above  pairwise  comparisons  of  a's  can  be 
made  by  the  linear  contrasts  of  a's.  The  probability  statement  for  the  family  of  all 
contrasts  is  given  by 


where 


k 

P  (|  Z  c • ( Y .  -  a. 
1  =  1  ' 

k 


<  c 


a  e 

k.k(n-l)  n 


k 

z 

i  =  l 


}  =  1  -  a. 


z  c-  =  0.  This  is  because  inequalities  (3.3)  hold  true  for  all 
i  =  l  1 


i^j  if  and  only 


if  the  inequality 


k  |c. 


|^ici(Yi  -  oi)|  1  ck?k(r!-l )  n  ^  2 


is  satisfied  for  all  linear  contrasts  I  c(  s  0.  Thus,  the  (l-a)100%  simultaneous  con- 

1  1 

fidence  intervals  for  all  linear  contrasts  of  the  location  parameters  is  given  by 
k  k  ~  k  |c.  | 

t“i“i  E  5  Ctn-l)  l  \z' 


These  comparisons  by  linear  contrasts  include  pairwise  comparisons  and  also 
the  comparisons  of  the  differences  between  weighted  averages  among  a's. 

Other  applications  of  the  C-di stribution  in  multiple  range  test  and  in  ranking 
and  selection  problem  could  be  expected. 
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4.  POWER  FUNCTION 


The  cumulative  distribution  function  (c.d.f. )  of  the  statistic  C  under  the 
alternative  hypothesis  :  a..  1  aj  may  be  derived  in  the  following  way: 

P(C  <  c| H1)  =  P(n(Y[)c]  -  Yfl])/e  <  c | Hj) 

k 

=  £  P(nY.  <  ce  +  nYj ,  y.  =  Y[ljt  for  all  i,  i  t  j | Hj > 


_  k 

£  P{mx(l,  +  6.., 
j=l  3  31 


0)  <  Z.  <  cX  +  1.  + 


for  all 


i  f  j). 


k  k 

=  £  r  r  H  {F(cx  +  Zi  +  6U)  -  FCmaxU, 
j=!°  °  i  =  l  J  31  3 


+  0))}  •  dF(2j)dG(x)  (4.1) 

where 

Zi  =  n(Yi  -  a..  )/9,  S-  =  n(aj  -  a.)/0  ,  for  i  =  1,  2,  k,  i^j, 

F(x)  =  1  -  e~x  and  6(x)  is  the  c.d.f.  of  X  =  e/e.  Note  that  for  any  given  i,  j,  the 
probability  F(cx  +  2j  +  '  F(max(zj  +  <5^,  0))  is  positive  if  cx  +  Zj  +  6-  >  0, 

zero  otherwise. 

The  c.d.f.  in  (4.1)  which  is  also  termed  the  probability  of  type  II  error  depends 
k-1 

on  the  (  2  )  differences  of  the  a's  and  it  cannot  be  simplified  further.  Thus,  the 
power  function  can  only  be  expressed  as  B(a)  =  1  -  p(C  <  c | ) ,  for  which  a  numerical 
integration  or  computer  simulation  should  be  used  under  It  is  clear  that  the  power 
function  under  riQ  reduces  to  (2.5)  as  a  special  case.  It  can  also  be  examined  that 


for  a  particular  configuration  of  ct's,  =  ...  *  a^,  ak  ■  +  6,  5 

(4.1)  can  be  reduced  to 

k*l 

P(c  <  c|H.)  =  e'(k‘1)d  L  (-l)£(k:1){(l  +  £c/v)~u 

1  e=o  1 

V  v-l  v.  ,  .  v-1-1  .  . 

"  {k_1)  TTvT  .?0  i!(  i  )(c}  (v+ic)-1'1 

z  ra!(J)  d’"  (k  - 

0=0  J  m=0  m  c 

where  v  =  k(n-l),  d  =  nfi/e,  c  = 

0  as  6'+oa. 


>  0,  expression 


,  which  converges  to  1-a  as  6-*0  and  to 


TABLL  1.  Upper  Percentage  Points  of  the  Range  Statistic  C 


2 

3 

4 

5 

6 

7 

8 

9 

10 

4.32 

4.81 

4.98 

5.08 

5.16 

5.23 

5.30 

5.36 

5.41 

3.11 

3.77 

4.09 

4.31 

4.48 

4.61 

4.73 

4.83 

4.92 

2.81 

3.48 

3.83 

4.08 

4.27 

4.42 

4.55 

4.66 

4.76 

2.67 

3.34 

3.71 

3.97 

4.16 

4.32 

4.46 

4.58 

4.68 

2.59 

3.26 

3.64 

3.90 

4.10 

4.27 

4.41 

4.53 

4.64 

2.54 

3.21 

3. 59 

3.86 

4.06 

4.23 

4.37 

4.50 

4.60 

2.50 

3.18 

3.56 

3.83 

4.04 

4.21 

4.35 

4.47 

4.58 

2.48 

3.15 

3.53 

3.81 

4.02 

4.19 

4.33 

4.46 

4.57 

2.46 

3.13 

3.52 

3.79 

4. 00 

4.17 

4.32 

4.44 

4.55 

2.43 

3.10 

3.49 

3.76 

3.98 

4.15 

4.  30 

4.42 

4.54 

2.39 

3.06 

3.46 

3.73 

3.95 

4.12 

4.27 

4.40 

4.51 

2.37 

3.04 

3.44 

3.71 

3.93 

4.11 

4.26 

4.39 

4.50 

2.35 

3.02 

3.41 

3.69 

3.91 

4.09 

4.24 

4.37 

4.48 

2.33 

2.99 

3.39 

3.67 

3.89 

4.07 

4.22 

4.35 

4.47 

2.31 

2.98 

3.38 

3.66 

3.88 

4.06 

4.21 

4.34 

4.46 

2.30 

2.97 

3.37 

3.65 

3.87 

4.05 

4.20 

4.34 

4.45 

18.00 

14.05 

12.23 

11.24 

10.64 

10.24 

9.96 

9.76 

9.61 

8.65 

8.43 

8.25 

8.14 

8.09 

8.06 

8.05 

8.06 

8.07 

6.93 

7.18 

7.27 

7.33 

7 . 39 

7.4  5 

7.51 

7.56 

7.62 

6.23 

6.64 

6.83 

6.96 

7.U7 

7.17 

7.25 

7.33 

7.40 

5.85 

6.34 

6.58 

6.76 

6.89 

7.01 

7.11 

7.19 

7.28 

5.61 

6.15 

6.43 

6.62 

6.77 

6.90 

7.01 

7.1] 

7.19 

5.45 

6.02 

6.32 

6.52 

6.69 

6.82 

6.94 

7.04 

7.13 

5.34 

5.92 

6.25 

6.45 

6.63 

6.77 

6.89 

7.00 

7.09 

5.25 

5.85 

6.17 

6.40 

6.58 

6.72 

6.85 

6. 96 

7.06 

5.12 

5.74 

6.08 

6.32 

6.51 

6.65 

6.79 

6.91 

7.01 

4.98 

5.62 

5.93 

6.23 

6.43 

6.59 

b  .73 

b .  85 

6.95 

4.90 

5.55 

5.92 

6.18 

6.38 

6.55 

6.69 

6.81 

6.9? 

4.79 

5 . 46 

5.84 

6 . 11 

6.32 

6.49 

6.64 

6.77 

6.88 

4.70 

5.38 

5.77 

6.05 

6.26 

6.44 

6 . 59 

6.72 

6 . 84 

4.65 

5.34 

5.73 

6.02 

6.24 

6.42 

b .  57 

6.70 

4.61 

5.30 

5.70 

5.99 

6.21 

6.39 

6.55 

6 . 68 

6.86 

TABLE  1*  ^Continued) 
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